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Introduction. In the context of the problem of plane deformation, a finite-element model of a natural landslide slope is 
developed. It allows for the joint work of a flow slide and a protective engineering structure. The Drucker-Prager model 
is used to take into account the physical nonlinearity of the slope layer material. To activate the kinematic instability, a 
viscoelastic interlayer is introduced into the design scheme, along which the landslide layer slides. 

Materials and Methods. Numerical experiments were performed using the ANSYS Mechanical software package, 
which implements the finite element method in the form of the displacement method. Slope discretization is performed 
on the basis of PLANE42 flat four-node finite elements. To simulate the displacement of the landslide layer relative to 
the fixed base, the combined viscoelastic elements COMBIN14 were used. 

Results. A physically nonlinear model of a natural landslide slope consisting of a base, a landslide layer, and a 
viscoelastic interlayer, 1s formalized. An engineering technique for analyzing the stress-strain state of the “slope- 
protective structure” system has been developed, taking into account the kinematic instability of the landslide layer. A 
series of computational experiments was carried out. 

Discussion and Conclusion. Based on the calculations performed, it is shown that the proposed method enables to 
specify the force action of the landslide layer on the protective structure and, thereby, to increase the reliability of the 


risk assessment when activating the landslide process. 
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Introduction. The design of anti-landslide engineering structures is one of the most important areas in the 
construction science. At the same time, numerical modeling of the stability of natural landslide slopes is associated with 
uncertainty factors, such as instability of the physical and mechanical characteristics of the material and the variability 
of the fracture trajectories. 

Nowadays, various variants of the limit equilibrium method and the method of reducing the shear strength of 
the material in combination with the finite element method are used to assess the slope stability [1-3]. The calculation 
result of the using the limit equilibrium method is the geometry of the “classical” circular sliding line. The method of 
reducing the strength provides more complete information about the sliding surface for the considered landslide area. 

Particularly, we can distinguish the direction of the study of slope processes based on the model of the 
centroidal motion of particles that elastically interact with each other and the inclined surface [4]. Another approach 
based on block kinematics is the so-called discontinuous deformation analysis (DDA) method [5]. The DDA method 


has limitations regarding the interaction between adjacent blocks. The essence of these restrictions is the impossibility 
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of mutual penetration of the blocks and the absence of tensile stresses between the blocks. Figure 1 (taken from [5]) 


shows the block coupling scheme used in the ADD method for modeling rockfalls. In this Figure, &,, k, and c,, c 


y 


— are the stiffness and viscosity parameters, respectively, in the direction of the axes _X and Y. 






Damper 


Fig. 1. Block diagram for modeling rockfalls 


The disadvantage of the DDA method is the impossibility of introducing engineering structures into the design 
scheme that prevent the progressive destruction of the mountain range. 

If, after assessing the natural slope stability, a decision is made on the feasibility of constructing landslide 
structures, then the data of engineering and geological surveys, as well as various scenarios of external force impacts on 
the landslide slope sections, are used as initial information for the design. Numerical simulation of the stress-strain state 
of a reinforced landslide slope is usually performed in a two-dimensional formulation within the framework of the 
nonlinear theory of elasticity using the finite element method. The ground pressure is assumed to be linearly distributed 
in depth. For engineering geotechnical calculations, specialized finite element software systems are used: Plaxis, 
MIDAS GTS, GEOS, which include the most common soil models. Modeling of complex geotechnical systems 1s 
carried out on the basis of the “heavy” software package ANSYS Mechanical [6]. 

It should be noted that in the papers devoted to the analysis of the stress-strain state of anti-landslide structures, 
there is no information on the use of models that take into account the direct sliding (creeping) of the landslide layer. In 
normal practice, design schemes are used, in which the joint (continuous) deformation of the entire soil mass 1s 
provided, up to the loss of the bearing capacity of the soil. This does not allow us to fully assess the force impact of the 
landslide layer on the structural elements. 

Materials and Methods. The design scheme of the landslide slope is shown in Fig. 2 (dimensions are given in 


meters). Consider that the slope consists of two layers: S, — base of the slope (dense clay); S, — landslide layer 


(waterlogged loam). We assume that the position of the sliding line is known. The external impact is presented by the 


self-weight of the landslide layer S', and the evenly distributed pressure gq . 





Fig. 2. Design scheme of landslide slope (dimensions are in meters) 


http://vestnik-donstu.ru 


The mechanical characteristics of the layer materials are given in Table 1. 
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Table 1 


Mechanical characteristics of slope layers 


Angle 
Deformation 


Poisson’s ratio of internal Density 


module £, 7 
friction 





Numerical simulation 1s performed using the ANSYS Mechanical software package. The finite element model 
of the slope is shown in Fig. 3. Two-dimensional 4-node finite elements PLANE42 are used for soil modeling. 


Moreover, the meshes of the layers S, and S', are topologically unrelated, 1.e., along the flow line, each node has a 


double numbering. This will ensure to model the possible displacement of the layer S, relative to the layer S 


2 LL? 

To describe the behavior of the slope layers, we use the Drucker-Prager model of an elastic-perfectly-plastic 
material [7], which is included in the ANSYS Mechanical complex. The equation of the yield surface for this model has 
the form 


F=1+3Bo,-o,=0, 

i m y 

where t, — intensity of tangential stresses; o,, — mean stress; B, o, — model parameters related to material 
constants C and 6 ratios: 


— 2sind _ 6Ccoso 


, 0, =. 
3(3-sin) (3(3- sino) 
Physical dependence linking stress and deformations for an elastic-perfectly-plastic material is described by the 


expression [8] 


Gy. 
oF a ; 
LV G./o. 
where y, — the intensity of shear deformations; G — the shear modulus; the tangential stress corresponding to the 


material tensile strength (Fig. 4) 
T,, =C+0, tego. 





Fig. 3. Finite element model of the slope Fig. 4. Graph T, ~ y, 


The effect of possible sliding of the layer S, relative to the layer S, is modeled using combined elements 


2 1 


(CE) of spring type COMBIN 14 (Fig. 5). 
The equilibrium equation for a combined spring element (CE) has the form 
[h]fiuj={p}, 
where /h/ — stiffness matrix; {u /, { p } — vectors of the columns of nodal displacements and forces. Depending 
on the orientation of the spring element, we have: 


— local CE axis € coincides with the axis x 
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1 0 -1 0 
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— local CE axis € coincides with the axis y 


0 0 0 0 

0 1 0 -1 fuj={0u,,0u,, } 5 
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0-10 1 PO =7 00.0 0..7 : 


Here, k,, k, — stiffness coefficients corresponding to the orientations of the spring along the axes x and y; 


I’ —symbol of the transpose operation. 





Fig. 5. Combined finite element COMBIN14 


We take the following values of the elasticity coefficients: 
k,=21 H/m, k,= 100 H/M. 


Viscosity parameters: c,= c,=0.5 N:s/m. 
When assigning coefficients &, and k,, the aim was to bring the model of the slope under consideration closer 
to reality. In particular, along the flow line, we consider the shear strain resistance and elastic rebound of layer S,.A 


similar approach to modeling the “slipping” effect was used in [9]. Note that in the accepted layer discretization, the 


nodes i and j, belonging to the spring elements coincide. The design load includes the self-weight of the layer S, 


and the evenly distributed force gq applied to the upper end of layer S', (Fig. 2). 


To organize the computational process, we use the nonlinear “solver” of the ANSYS Mechanical complex, 
which implements the Newton-Raphson method. 

To verify the finite element model, we estimate the slope stability under consideration through excluding the 
COMBIN14 elements from the design scheme and thus obtaining an unreduced model. By calculation, it is established 


that for an unweakened slope, the pressure limit value q , for the given mechanical characteristics of layers S, and S,, 


is 400 kPa. The operation process of the nonlinear solver is shown in the form of a graph in Fig. 6. Here, along the 
abscissa axis, the total number of iterations at each step of “Time” loading is plotted, along the ordinate axis, the 
corresponding norms of the residual forces “F” and displacements “U” are plotted. The “Time” parameter ranges from 0 
to 1. The iterative process continues until the values “F’—“L2” are less than the values “F’—“CRIT”’. The intersection of 
the peak of the sawtooth curve “F’”—“L2” and the curve “F’”—“CRIT” indicates that, in this example, the iterative process 
“converged” with the specified accuracy. 


Figure 7 shows a picture of the distribution of horizontal displacements u, in the body of an unweakened 


slope at g =400 kPa. As follows from the Figure, the loss of stability of the slope is local in nature. In this case, the 
circular-cylindrical type of the flow line is clearly traced. The corresponding distribution of the intensity field of plastic 


shear deformation y, is shown in Fig. 8. 
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Fig. 6. Convergence graph of the Newton-Raphson procedure 
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Fig. 7. Visualization of displacements u, at g =400 kPa 





Fig. 8. Visualization of distribution y, at g =400 kPa 


Research Results. To simulate the effect of kinematic instability of the landslide layer S,, we analyze the 
stress-strain state of the slope taking into account the weakening at different pressure values g. We emphasize that in 


this case, the elements of COMBIN]4 are included. Fig. 9 shows the graph u,~gq, where u, — the horizontal 


displacement of point J of the “tongue” of the landslide layer S, (Fig. 2). 


Comparing the results of calculations for weakened and unweakened slopes, we come to the conclusion that 


the inclusion of a sliding layer in the design model significantly (by more than one order of magnitude) reduces the 


limit value q. 
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Fig. 9. Graph of the dependence of horizontal movement of point J of the “tongue” of landslide layer S, on pressure value q 


The graph uw. ~gq _ shows that, for the considered slope model, the landslide process has a pronounced 


nonlinear character. Moreover, in the pressure range from 6 kPa to 11 kPa, there is a displacement of the “tongue” 


of layer S, in the opposite direction to the expected landslide process. Fig. 10 shows a contrasting pattern of the 


2 


distribution of the displacement field in the x- axis direction corresponding to the loading of layer S$, with 


pressure g =16.5 kPa. 
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Fig. 10. Visualization of displacements uw, at gq =16.5 kPa 


A fragment of a finite element model of a landslide “tongue” in a deformed state at g =16.5 kPa is shown in 


Fig. 11. 





=, 


Fig. 11. Visualization of a finite element model of a slope fragment in a deformed state at g =16.5 kPa 


As can be seen from Fig. 11, the proposed method allows you to simulate the effect of sliding of layer S, 
relative to layer S,. This effect is shown in the relative displacement of the grid nodes (shown by the arrow). Maximum 


relative displacement u reaches 18.2 cm, which, from a physical point of view, 1s quite consistent with the landslide 


xX max 


process onset. 


Fig. 12 shows the considered variants of schemes of anti-landslide structures made of concrete. The following 


mechanical constants for concrete are accepted: E.= 3-10* MPa; vV,= 0.2; p, = 2446 kg/m’. The design value of the 


height of the protective structures 4 =2.6 m. The depth of the pile foundation is 8 m (schemes 5 and c). Pressure value 
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g = 16.5 kPa. 


The pile foot (scheme b) and the sheet piling (scheme c) are presented by 2-node beam elements BEAM3. The 
diameter of the piles is assumed to be 0.35 m. 
The visualization of the considered schemes of anti-landslide structures in the deformed state is shown in 


Fig. 13. From the above data, it follows that when modeling the landslide process, the “tongue” of layer S', 1s pressed 


2 


into the foundation S, This explains the “collapse” of anti-landslide structures towards the movement of layer S 


h Yo h 
sa $e 
a — 

h Yo h 
ad $e 


Fig. 12. Schemes of anti-landslide structures: 


2° 






a) retaining wall, b) retaining wall on a pile foundation, c) sheet piling, d) gabion 


Graphs of changes in horizontal displacements wu, and stresses o, along the height A of the considered 


protective structures are shown in Fig. 14. 

Analyzing the results presented, we come to the conclusion that the most pronounced anti-landslide effect from 
an engineering point of view is scheme b, which provides a relatively uniform horizontal displacement of the retaining 
wall without “collapse’’. In this case, maximum horizontal displacement of the wall is 16 mm. To increase the reliability 
of this design, you can use anchor ties [10]. 

A certain practical interest in the design of anti-landslide structures is the nature of the distribution of internal 
forces in the pile foundation. For scheme 5, the diagrams of the bending moments M occurring in the left and right 


piles are shown in Fig. 15. The value of the maximum normal compression stress under bending for scheme b was o 


ax 


= 337 kN/m’, which is significantly less than the compressive strength of concrete (25.5 MPa). 
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c) d) 


Fig. 13. Visualization of protective structures in the deformed state: 


a) retaining wall, b) retaining wall on a pile foundation, c) sheet piling, d) gabion 
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Fig. 14. Graphs u, ~ A and o, ~ A for protective structures: 


a) retaining wall, b) retaining wall on a pile foundation, c) sheet piling, d) gabion 
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Fig. 15. Diagrams of bending moments in piles of scheme 6 in Fig. 12: 


a) left pile; 5) right pile 
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Diagrams of shear forces Q in piles are shown in Fig. 16. 
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Fig. 16. Diagrams of shear forces in piles of scheme 6 in Fig. 12: 


a) left pile; 5) right pile 


Fig. 16 shows that the cross-section of a possible pile cut is located below the zero mark of the structure. 

Discussion and Conclusions. Under the problem of plane deformation, a finite element model is constructed 
that simulates the landslide process of a natural slope. In contrast to the existing methods of calculating landslide 
structures, this approach implements the concept of kinematic instability of the landslide layer, the essence of which is 
to use combined finite elements along the assumed sliding line. A series of numerical experiments was used to test the 
proposed method for studying the interaction of a sliding landslide and various types of protective structures. It is 
established that the most effective anti-landslide protection at the given geometric and physical parameters 1s provided 


by a combined structure consisting of a retaining wall and a pile foundation. 
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